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Abstract 

The electronic stress tensor of the hydrogen molecule ion H^" is investigated for the ground state 
(<7 9 ls) and the first excited state (<r*ls) using their exact wave functions. A map of its largest 
eigenvalue and corresponding eigenvector is shown to be closely related to the nature of chemical 
bonding. For the ground state, we also show the spatial distribution of interaction energy density 
to describe in which part of the molecule stabilization and destabilization take place. 
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The stress tensors are used widely for description of internal forces of matter. They have 
been also established in quantum mechanical context pQ, although it seems that basic idea 
dates back to Ref. [2]. Later, applications for quantum systems are investigated in many 
literatures from various aspects [31 IH El El El El El HDl HH H21 LTHl HH UHl HEl HZl HE! UHl 
1201 [211 122] . For example, Ref. [6] and followers have focused on stress tensor, which is 
associated with forces on nucleus. Our focus in this paper is electronic stress tensors, which 
measure effects caused by internal forces acting on electrons in molecules and particularly 
those between bonded atoms [121 HSl US]- Such electronic stress tensor has been successfully 
used to define bond orders expressing the strength of chemical bonds [18J. It should be 
noted that people sometimes use different definitions for the stress tensor due to their own 
purpose and formalism. We use the expression derived by one of the authors in Ref. [12] 
which we believe most straightforward and suitable to discuss chemical bonds. 

In this paper, we elucidate the role of the electronic stress tensor using hydrogen molecule 
ion H^~, the simplest of all molecules, whose exact wave functions are known. Although 
there have been analyses of electronic stress tensor of other molecules using approximate 
wave functions expanded in the Gaussian functions [121 [El [161 LTHl EH G2] , analysis using 
exact wave functions has not been conducted. We calculate spatial distributions of stress 
tensor and other associated quantities from the exact wave functions of and examine 
their detailed features that can be used to illustrate the chemical bond of the molecule. We 
believe that this paper demonstrates usefulness of the electronic stress tensor for an analysis 
of chemical bonding in the simplest but most exact form. 1 

We begin by showing our expression of the electronic stress tensor rfj constructed from 
the wave function ip. That is, 

4 = ^ W - StfW + cc.) , (1) 

where {i, j} = {1, 2, 3} denote spatial coordinates, m is the electron mass, and cc. stands for 
complex conjugate. The expression for many-electron systems, that is, for almost every other 
molecule, is found in Ref. [12]. There, the stress tensor has been derived by field theoretic 
method which is applicable to many-particle system. Here, for illustrative purpose, we derive 
it using the Schrodinger equation for a single electron in H^". 

1 Calculations of stress tensor field as regards a H^" molecule are found in Refs. [23] but with different 
definitions from ours. 
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The time-dependent Schrodinger equation reads 

where the potential energy is assumed to be real. It is well-known that the continuity 
equation dn/dt + V • j — holds by defining probability density n = \ip\ 2 and probability 
flux j = (h/m)lm(il)*'Vip). Now, the equilibrium equation for electronic stress tensor is 
obtained from an equation of motion of j. Namely, using Eq. (pj, time derivative of j can 
be written as 

% = --^ (W ■ W - fVVV + c.c.) - ^ W. 
at 4m z v ' m 

For the steady state, which is the case of molecular systems we consider, since dj/dt = 0, 
we obtain the equilibrium equation 

djrfj + F Lti = 0, (3) 

where Fn = — \ip\ 2 diV is the Lorentz force generated by nuclei. For later convenience, we 
here define the first term in the equation as tension F Tti = djrf-. From this equation, we see 



J 3 •ij- 

that rfj is the stress caused by purely quantum mechanical effect. 



We compute the stress tensor defined by Eq. ([T]) using the exact wave functions of H^" [25] . 
As is usually done for the stress tensor, we examine its largest eigenvalue and corresponding 
eigenvector. The sign of the largest eigenvalue tells whether electrons at a certain point in 
space feel tensile force (positive eigenvalue) or compressive force (negative eigenvalue) and 
the eigenvector tells direction of the force. 

The case for the ground state (er s ls) is plotted in the plane including two H nuclei as 
shown in the left panel of Fig. [l] The origin of the coordinate is taken to be the midpoint 
of the two nuclei. The eigenvalue along the internuclear axis is also plotted as red solid line 
in Fig. [2j The internuclear distance is taken to be 2.0 bohrs, which gives the equilibrium 
distance R e (to be more precise, R e = 1.9972 bohrs j2S] but it does not make practical 
difference in our argument below). The region with positive eigenvalue spreads between 
the nuclei which corresponds to the formation of a chemical bond. The electronic stress 
tensor has this useful feature that it can express something is pulled up (that is, tensile 
force exists) when a covalent bond is formed. In detail, the positive eigenvalue region is 
bounded by a closed sphere-like surface that touches two H nuclei. Also, in that region, the 
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FIG. 1: The spatial distributions of the largest eigenvalue and corresponding eigenvectors of stress 
tensor of molecule are plotted in the plane including two H nuclei. They are located at 
(z,x) = (—1.0,0.0) and (1.0, 0.0). The left panel is for the ground state (<r ff ls) and the right panel 
is for the first excited state (a* Is). The black solid line shows the zero surface of the eigenvalue. 
For dgls state, it is positive inside the circle and is negative outside. For the <r„ls state, it is 
negative on both sides of the surface. 
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FIG. 2: The largest eigenvalue of the stress tensor along the internuclear axis (along the z-axis in 
Fig. [T]). The red solid line is for the ground state (cr g ls) and the green dashed line is for the first 
excited state (<r*ls). 
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FIG. 3: The spatial distributions of the tension F T of molecule are plotted similarly to Fig. |l| 
The left panel is for the ground state (<r ff ls) and the right panel is for the first excited state (<r*ls). 
Normalized tension vectors F T /\F T \ are shown by arrows and the norm \F T \ is depicted by a color 
map. |-F T | diverges at the positions of the nuclei. Contours which correspond to \F T \ = 5 are 
drawn by black dashed lines. \F T \ = at (0.0, 0.0) for the <r 9 ls state and on z = for the cr*ls 
state, where arrows are not shown. 

eigenvectors forms a bundle of flow lines that connects the H nuclei. Such a region, called 
"spindle structure" [15], is clearly visible using the analysis with the exact wave function. 

The case for the first excited state (<r*ls) is shown in the right panel of Fig. [I] on the 
same plane and with the same internuclear distance. See also the green dashed line in Fig. [2] 
for the eigenvalue along the internuclear axis. Since this state does not form a bound state, 
we do not expect a bond between the H nuclei. Indeed, there is no region with positive 
eigenvalue between the H nuclei for this case. In addition, eigenvectors are perpendicular to 
the internuclear axis at the region around the middle of the two H nuclei, showing there is 
no force connecting them. 

These results, positive eigenvalue for bonding orbital and negative eigenvalue for anti- 
bonding orbital, are consistent with Ref. [15] . where these features are pointed out using H 2 
and He2, respectively. 

We now examine the tension field F T ^ = djrfj of the H^" molecule. Writing explicitly, 

h 2 

F T = (W • VV - V*VV 2 V + c.c.) . (4) 

4m v ' 
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Since this force, which cancels the classical Lorentz force at each point in space [Eq. ([3j], 
also expresses purely quantum mechanical effects, it is considered to carry some information 
of a chemical bond. The spatial distribution of tension vector is plotted for <J g ls and ex* Is 
states in Fig. [3j On comparing Figs. |1] and [3} we see that in contrast with the case of the 
eigenvalue and eigenvector of the stress tensor, there is not much difference between the 
a g ls state and the a* Is state for the case of the tension vector. In Fig. [TJ there is a striking 
difference between cx ff ls and ex* Is that the former develops a region with positive eigenvalue 
while the latter exhibit negative region for the most of the space. There is also a marked 
difference in the pattern of eigenvectors at the internuclear region. In Fig. |3j the two states 
show a relatively similar pattern of the tension vector field. One difference can be pointed 
out that there is a point (x, z) = (0.0, 0.0) at which tension vanishes for the cr g ls state while 
the a* Is states has a plane z = on which it vanishes. However, overall feature resembles 
closely to each other. This indicates that the stress tensor has much more information as 
regards a chemical bond. 

We note, however, that the tension can provide a handy way to characterize a chemical 
bond. Actually, it has been used to define an index of bond strength and turned out to 
be very useful [TH [191 E2]. In Ref. [H], it has been proposed to calculate energy density 
[defined below, Eq. ([6])] at the "Lagrange point," which is defined as the vanishing point of 
the tension such as we have found at the midpoint of the internuclear axis of the ground 
state H^" molecule in Fig. [3] It should also be noted that the existence of the Lagrange 
point alone cannot guarantee a chemical bond as is easily seen from the tension of the ex* Is 
state of Hj, which has the Lagrange point (or rather the Lagrange plane z — 0). Hence, the 
tension can be a useful tool when combined with analysis of the stress tensor. 

Finally, we describe our definition of energy density, which is derived from the stress 
tensor, and discuss the interaction energy density of the ground state of the H^" molecule. 
In Ref. [12J, it has been proposed to define energy density e T from the trace of the stress 
tensor as 

er = \Y,4 (5) 

i 

h 2 

= — U*V 2 ij - • + c.c.) . (6) 

8m v ' 

Note that this definition gives correct total energy when integrated over the whole space 
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FIG. 4: The spatial distributions of the interaction energy density Ae r (left panel) and the dif- 
ferential electron density An (right panel) of H^" molecule are plotted for the ground state (<r 9 ls). 
The zero surface is shown by the red solid line for each panel. The interaction energy diverges at 
the positions of the nuclei and contours that correspond to e T = —0.1 are shown by white dashed 
lines in the left panel. 

and the virial theorem is applied. Furthermore, we consider the difference of the energy 
density between the ground state and the separated atoms. We call this interaction energy 
density Ae T [18], and it is a useful quantity to discuss where in a molecule stabilization or 
destabilization take place by chemical bond formation. 
For the case of H^", it is appropriate to define Ae T as 

Ae T = £ r (rFJ) , (7) 

where Ha denotes the H atom at z = —1.0 and Hb at z = 1.0 (as noted above, more precise 
R e is 0.3% shorter and the virial theorem holds only at the equilibrium position but we can 
practically consider R e to be 2.0 bohrs). We plot this in the left panel of Fig. |4j From the 
definition, negative Ae T (shown in blue) indicates stabilized regions and positive Ae T (shown 
in red) indicates destabilized regions. We see that stabilization mostly takes place between 
the H nuclei. The stabilized region slightly expands outside of the internuclear region and 
the zero surface of Ae T is found to include the two nuclei. 
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It is instructive to compare Ae T with the differential electron density An, 

Au = n { nt)- n ^ )+ 2 n ^\ (8) 



where n = \i/j\ 2 is the ordinary electron density. This is plotted in the right panel of Fig. |4j 
We see that there is a region with increased electron density (shown in red) between the 
H nuclei and one with decreased electron density (shown in blue) outside the internuclear 
region. Compared with the map of Ae T , the region with increased electron density corre- 
sponds nicely to the stabilized region, which is consistent with a conventional picture of 
chemical bond formation. However, it should also be noticed that Ae T and An do not have 
one-to-one correspondence. For example, we see in Fig. [4] that Ae T has slightly larger zero 
surface that that of An. Moreover, integrating An over the whole space gives zero whereas 
Ae T gives the difference of binding energy between and H atom. 

In conclusion, we have calculated the stress tensor of the simplest molecule, H^~, and 
demonstrated its relevance to chemical bonding. It is not only the simplest but is also 
most accurate due to the existence of the exact wave functions. By using them, we are 
able to investigate the eigenvalue and eigenvector of the stress tensor and the tension vector 
field of both the ground state and the first excited state in a very accurate manner. It 
also makes possible the precise comparison between the interaction energy density, which 
tells the stabilized/destabilized regions, and conventional differential electron density with 
respect to the ground state. We believe that this calculation can be very useful as the most 
basic example which shows an important role of the electronic stress tensor for illustrating 
a chemical bond. Moreover, its exactness would serve as a benchmark to test validity of the 
calculation of the stress tensor which uses approximate wave functions. This is also very 
important since we can only obtain approximate wave functions for almost every molecular 
system. Therefore, future investigation should include how the stress tensor of Hj" deviates 
from the exact form as we change the level of approximation of the quantum chemistry 
calculation to derive H^" wave functions. Such a study will solidify the basis of the stress 
tensor analysis of molecular systems and open up areas of further application. 
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